Method and device for extracting skeleton topology structure of electric power grid

ABSTRACT

A system and method for extracting a skeleton topology structure for an electric power grid, the method comprising: receiving a description of a topology sub-structure corresponding with user&#39;s need and a description of skeleton topology sub-structure extracted from the topology sub-structure; generating a first incidence matrix based on the description of the topology sub-structure and a second incidence matrix based on the description of the skeleton topology sub-structure; generating a third incidence matrix based on a primary topology structure of electric power grid; searching from the third incidence matrix a sub-matrix that matches the first incidence matrix; obtaining a fourth incidence matrix by using the second incidence matrix to transform the matching sub-matrix; and generating a skeleton topology structure corresponding to the primary topology structure based on the fourth incidence matrix.

CROSS-REFERENCE TO RELATED APPLICATIONS

The present invention relates to and claims the benefit of the filingdate of commonly-owned, PCT Patent Application No. PCT/CN2012/074087,filed Apr. 16, 2012, which further claims the benefit of priority dateof commonly-owned, co-pending Chinese Patent Application No. CN201110140427.8, filed on May 27, 2011, the entire contents anddisclosure of which is incorporated by reference as if fully set forthherein.

TECHNICAL FIELD

The present invention generally relates to an electric power grid, andmore specifically to a method and device for extracting skeletontopology structure of the electric power grid.

DESCRIPTION OF THE RELATED ART

An electric power system as a whole consists of electric componentsconnected in a certain form, such as generators, transformers, buses,switches, disconnector, and lines. Performance of its electric functionis restricted by component characteristics and connecting relations ofthe components. With regardless of physical component characteristics ofthe grid, the electric power system can be abstracted as a grid topologyformed by branches and nodes connected by the branches. Considering therestrictions of component characteristics and their connectingrelations, an electric power grid actually includes two types oftopology structures: geometry topology and physical topology. Thegeometry topology represents the geometrical connecting status of griddevices, while the physical topology represents the physical electriccoupling relations of grid components.

Analysis of grid topology structure of an electric power system isusually conducted in the following two aspects. The first aspect lies inestablishing of the grid topology structure. According to the status ofswitches, the grid connected by various devices can be indicated by amatrix of connecting relations of branches and nodes, which can be usedfor analyzing and computing the power system. Also, identification ofinter-connected individual sub-device is the ground for physicalanalyzing, computing and studying of the electric power system. Thesecond aspect lies in studying and utilizing the grid topology structureto excavating the internal relations between its topology structure andphysical functions, so as to make it more convenient and easier foranalyzing and controlling the electric power system. In actualoperation, different business departments of an electric power gridcompany need to make analysis and computing of the electric power systemaccording to their own business needs. Thereby, different points ofconcerns to primary topology structure of the electric power system arefocused on abstracting different skeleton topology networks from aprimary topology structure. Currently, the business departments of theelectric power grid company mainly extracts skeleton topology structurefrom the primary topology structure manually according to their businessneeds and based on experience of experts. However, there are actually ahuge number of components within an electric power system and the linesare heavy and complicated. For reasons such as load changes, equipmentmaintenance and active optimizations, the grid is continuously renovatedat the equipment level, and the primary topology structure is oftenchanged. If the business departments extract and maintain skeletontopology structure from the primary topology structure manually, theworking procedures would be overloaded with details and the costs wouldbe high.

Therefore, it is necessary to provide a method for automatic andefficient extracting skeleton topology structure according tooperational needs.

SUMMARY OF THE INVENTION

In order to solve these problems, this invention provides an automaticand efficient method and device for extracting skeleton topologystructure.

According to one aspect of the invention, a method for extracting askeleton topology structure of an electric power grid is provided. Themethod comprises: receiving a description of a topology sub-structurecorresponding with user's need and a description of a skeleton topologysub-structure extracted from the topology sub-structure; generating afirst incidence matrix based on the description of the topologysub-structure and generating a second incidence matrix based on thedescription of the skeleton topology sub-structure; generating a thirdincidence matrix based on a primary topology structure of electric powergrid; searching from the third incidence matrix a sub-matrix thatmatches the first incidence matrix; obtaining a fourth incidence matrixby using the second incidence matrix to transform the matchingsub-matrix; and generating a skeleton topology structure correspondingto the primary topology structure based on the fourth incidence matrix.

According to a second aspect of the invention, a device for extracting askeleton topology structure of an electric power grid is provided. Thedevice comprises: a receiving module configured to receive a descriptionof a topology sub-structure corresponding with user's need and adescription of a skeleton topology sub-structure extracted from thetopology sub-structure; an incidence matrix generating module configuredto generate a first incidence matrix based on the description of thetopology sub-structure, to generate a second incidence matrix based onthe description of the skeleton topology sub-structure, and to generatea third incidence matrix based on a primary topology structure ofelectric power grid; a searching module configured to search from thethird incidence matrix a sub-matrix that matches the first incidencematrix; a matrix transforming module configured to obtain a fourthincidence matrix by using the second incidence matrix to transform thematching sub-matrix; and a skeleton topology structure generating moduleconfigured to generate a skeleton topology structure corresponding tothe primary topology structure according to the fourth incidence matrix.

The method and device for extracting skeleton topology structureaccording to this invention can automatically and efficiently extractskeleton topology structure according to business needs.

BRIEF DESCRIPTION OF THE DRAWINGS

In combination with the drawings and in reference to the followingdetailed descriptions of exemplified embodiments, the invention itself,preferred embodiments, objectives and advantages of this invention canbe better understood.

FIG. 1 shows the method for extracting a skeleton topology of theelectric power grid according to an embodiment of this invention.

FIG. 2 shows a diagram of the topology substructure, according to anembodiment of this invention.

FIG. 3 shows a directed graph of the connecting relations ofnodes-branches formed by abstracting the topology substructure in FIG.2, according to an embodiment of this invention.

FIG. 4 shows diagram of the skeleton topology substructure, according toan embodiment of this invention.

FIG. 5 shows a directed graph of the connecting relations ofnodes-branches formed by abstracting the skeleton topology substructurein FIG. 4, according to an embodiment of this invention.

FIG. 6 shows diagram of the primary topology structure, according to anembodiment of this invention.

FIG. 7 shows a directed graph of the connecting relations ofnodes-branches formed by abstracting the topology substructure in FIG.6, according to an embodiment of this invention.

FIG. 8 shows a flow chart of searching an incidence matrix of thetopology substructure within the incidence matrix of a primary topologystructure of the electric power grid.

FIG. 9 shows a specific process for implementing Step S802 in FIG. 8.

FIG. 10 shows a specific process for implementing Step S804 in FIG. 8.

FIG. 11 shows a skeleton topology structure corresponding to the primarytopology structure, according to an embodiment of this invention.

FIG. 12 shows a block diagram of the device for extracting skeletontopology in an electric power grid, according to an embodiment of thisinvention.

FIG. 13 is a schematic block diagram of the structure of the computingdevices capable of implementing the embodiments of this invention.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT

In the following discussion, the method and device for extracting askeleton topology structure of an electric power grid according to theembodiments of this invention will be described in combination with theattached drawings, so that objectives and advantages of this inventioncan be better understood.

FIG. 1 shows a method for extracting a skeleton topology of the electricpower grid according to an embodiment of this invention. The methodcomprises: in Step S101, receiving descriptions of a topologysub-structure corresponding with user's need and receiving descriptionsof skeleton topology sub-structure extracted from the topologysub-structure; in Step S102, generating a first incidence matrix basedon descriptions of the topology sub-structure and generating a secondincidence matrix based on descriptions of the skeleton topologysub-structure; in Step S103, generating a third incidence matrix basedon a primary topology structure of the electric power grid; in StepS104, searching from the third incidence matrix a sub-matrix thatmatches the first incidence matrix; in Step S105, generating a fourthincidence matrix by using the second incidence matrix to transform thematching sub-matrix; in Step S106, generating skeleton topologystructure corresponding to the primary topology structure according tothe fourth incidence matrix.

In step S101, a description of a topology sub-structure that correspondswith the need of user and a description of a skeleton topologysub-structure extracted from the topology sub-structure are received.Because different business departments of an electric power company havedifferent business needs and different business analysis objectives,their points of concerns are different even for a same primary topologystructure of electric power grid. Thereby, different businessdepartments need to specify topology substructures that correspond withtheir business needs, and specify skeleton topology substructuresextracted from the topology substructures. With regards to a primarytopology structure of a whole electric power system, the topologysubstructure is subset of the primary topology structure, and theskeleton topology sub-structure is further simplification of thetopology substructure, which is made by the departments according totheir business needs. According to an embodiment of this invention, theuser is not only provided with an editing tool with graphic interface sothat the user can edit topology structures concerned personally, butalso can select a topology substructure concerned from a primarytopology structure and make further simplification to the topologysubstructure to form a skeleton topology substructure. According toanother embodiment of this invention, files containing topologysubstructures and skeleton topology substructures edited by the user canbe imported.

In step S102, a first incidence matrix is generated based on thedescription of the topology sub-structure, and a second incidence matrixis generated based on the description of the skeleton topologysub-structure, in which the first incidence matrix corresponds to theincidence matrix of the topology sub-structure and the second incidencematrix corresponds to the incidence matrix of the skeleton topologysub-structure. In the field of electric power transmission, theincidence matrix can be used to describe an incidence relationshipbetween the connections among electric components and the electriccomponents within a topology structure of the electric power grid.Specifically, the primary topology structure of an electric power gridconnecting various devices can be abstracted as a directed graph of theconnecting relations of nodes-branches, in which the nodes indicate theelectric components while the branches indicate the connections of theelectric components. A topology structure of an electric power grid withn nodes and b branches can be indicated by an n×b matrix. Row vectors ofthe matrix correspond to the nodes while column vectors of the matrixcorrespond to the branches of the electric power grid. Vector a_(ij) thematrix is defined as:

$a_{ij} = \left\{ \begin{matrix}{{0\left( {{when}\mspace{14mu}{branch}\mspace{14mu} k\mspace{14mu}{is}\mspace{14mu}{not}\mspace{14mu}{connected}\mspace{14mu}{to}\mspace{14mu}{node}\mspace{14mu} j} \right)};} \\{{1\begin{pmatrix}{{{when}\mspace{14mu}{branch}\mspace{14mu} k\mspace{14mu}{is}\mspace{14mu}{connected}\mspace{14mu}{to}\mspace{14mu}{node}\mspace{14mu} j}\mspace{11mu}} \\{\;{{in}\mspace{14mu}{the}\mspace{14mu}{direction}\mspace{14mu}{that}\mspace{14mu}{the}\mspace{14mu}{current}\mspace{14mu}{flows}\mspace{14mu}{out}\mspace{14mu}{from}\mspace{14mu}{node}\mspace{14mu} j}}\end{pmatrix}};} \\{{- 1}\begin{pmatrix}{{{when}\mspace{14mu}{branch}\mspace{14mu} k\mspace{14mu}{is}\mspace{14mu}{connected}\mspace{14mu}{to}\mspace{14mu}{node}\mspace{14mu} j}\mspace{14mu}} \\{{in}\mspace{14mu}{the}\mspace{14mu}{direction}\mspace{14mu}{that}\mspace{14mu}{the}\mspace{14mu}{current}\mspace{14mu}{flows}\mspace{14mu}{into}\mspace{14mu}{node}\mspace{14mu} j}\end{pmatrix}}\end{matrix} \right.$

Each column of the incidence matrix corresponds to a branch, everybranch must connect two nodes, and the directions are one inward and theother outward. Thus, every column contains only two nonzero elements, +1and −1. If summing elements of all rows by columns, a row with eachelement as zero would be obtained. Thereby, the rows of a matrix are notindependent to each other. For elements of any of the rows, they can beobtained by summing all rows other than this row and changing the signsof the sum. Taking the topology substructure in FIG. 2 as an example todescribe the generation of corresponding incidence matrix, the firststep is to abstract the topology substructure in FIG. 2 to form thedirected graph of the connecting relations of nodes-branches. It can beseen from the directed graph that the connection between node 1 and node3 is disconnector 1, the connection between node 2 and node 4 isdisconnector 2, the connection between node 3 and node 5 is grounddisconnector 1, the connection between node 4 and node 5 is grounddisconnector 2, the connection between node 3 and node 4 is breaker 1. Arow vector V includes bus segment 1, bus segment 2, dummy point 3, dummypoint 4, dummy point 5. A column vector M includes disconnector 1,disconnector 2, breaker 1, ground disconnector 1, ground disconnector 2.According to the above description of the incidence matrix, an incidencematrix A=V×M corresponding to the topology substructure is obtained.

${{Row}\mspace{14mu}{vector}\mspace{14mu} V} = \begin{bmatrix}{{Bus}\mspace{14mu}{segment}\mspace{14mu} 1} \\{{Bus}\mspace{14mu}{segment}\mspace{14mu} 2} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 3} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 4} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 5}\end{bmatrix}_{5*1}$${{Column}\mspace{14mu}{vector}\mspace{14mu} M} = \begin{bmatrix}{{Disconnector}\mspace{14mu} 1} \\{{Disconnector}\mspace{14mu} 2} \\{{Breaker}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}_{5*1}^{T}$${{Incidence}\mspace{14mu}{matrix}\mspace{14mu} A} = \begin{bmatrix}1 & 0 & 0 & 0 & 0 \\0 & 1 & 0 & 0 & 0 \\{- 1} & 0 & 1 & 1 & 0 \\0 & {- 1} & {- 1} & 0 & 1 \\0 & 0 & 0 & {- 1} & {- 1}\end{bmatrix}_{5*5}$

Again taking the topology substructure in FIG. 2 as an example, becausethe user is a dispatch department of electric power system, moreconcerns are given to the status of buses and switches, the skeletontopology substructure extracted according to user needs only includebuses and switches. FIG. 4 shows the skeleton topology substructureaccording to an embodiment of this invention. FIG. 5 shows a directedgraph of the connecting relations of nodes-branches formed byabstracting the skeleton topology substructure in FIG. 4, according toan embodiment of this invention. With the method of generating theincidence matrix, an incidence matrix A′ corresponding to the skeletontopology substructure is obtained.

${{Incidence}\mspace{14mu}{matrix}\mspace{14mu} A^{\prime}} = \begin{bmatrix}0 & 0 & 1 & 0 & 0 \\0 & 0 & {- 1} & 0 & 0 \\0 & 0 & 0 & 0 & 0 \\0 & 0 & 0 & 0 & 0 \\0 & 0 & 0 & 0 & 0\end{bmatrix}_{5*5}$

In step S103, a third incidence matrix is generated based on a primarytopology structure of the electric power grid, in which the thirdincidence matrix corresponds to the incidence matrix of the primarytopology structure. The primary topology structure of the electric powergrid corresponds to the electric power grid as a target of analysis.Generally, the electric power grid, as a target of analysis, covers acity or a district. A city-level medium electric transmission system hasabout 400 lines, 120 substations. Each substation has about 500 primarydevices and each line has about 50 primary devices. Major devicesinclude: main transformers, switches, breakers, buses, lines, currenttransformers, voltage transformers, reactors and ground disconnector.That means the primary topology structure of corresponding electricpower grid has complicated lines and numerous electric components. Basedon the above description of the incidence matrix, the incidence matrixof primary topology structure is generated according to the primarytopology structure of the electric power grid. To make it convenient fordescriptions, a simple primary topology structure is schematicallyselected in an embodiment of this invention. FIG. 6 shows a diagram ofthe primary topology structure, according to an embodiment of thisinvention. By abstracting the topology substructure in FIG. 6, thedirected graph of the connecting relations of nodes-branches in FIG. 7is formed. The row vector V′ includes auxiliary bus segment 1, auxiliarybus segment 2, dummy point 3, dummy point 4, dummy point 5, main bussegment 2, main bus segment 1, outgoing line 2, outgoing line 1. Thecolumn vector M′ includes disconnector 1, disconnector 2, breaker 1,ground disconnector 1, ground disconnector 2, breaker 2, breaker 3,disconnector 3, disconnector 4. Then according to the description of theincidence matrix, the incidence matrix A″=V′×M′ corresponding to theprimary topology structure is obtained.

$\mspace{20mu}{{{Row}\mspace{14mu}{vector}\mspace{14mu} V^{\prime}} = \begin{bmatrix}{{Auxiliary}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 1} \\{{Auxiliary}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 2} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 3} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 4} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 5} \\{{Main}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 2} \\{{Main}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 1} \\{{Outgoing}\mspace{14mu}{line}\mspace{14mu} 2} \\{{Outgoing}\mspace{14mu}{line}\mspace{14mu} 1}\end{bmatrix}_{9*1}}$$\mspace{20mu}{{{Column}\mspace{14mu}{vector}\mspace{14mu} M^{\prime}} = \begin{bmatrix}{{Disconnector}\mspace{14mu} 1} \\{{Disconnector}\mspace{14mu} 2} \\{{Breaker}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2} \\{{Breaker}\mspace{14mu} 2} \\{{Breaker}\mspace{14mu} 3} \\{{Disconnector}\mspace{14mu} 3} \\{{Disconnector}\mspace{14mu} 4}\end{bmatrix}_{9*1}^{T}}$${{Incidence}\mspace{14mu}{matrix}\mspace{14mu} A^{''}} = \begin{bmatrix}1 & 0 & 0 & 0 & 0 & 1 & 0 & 0 & 0 \\0 & 1 & 0 & 0 & 0 & 0 & 1 & 0 & 0 \\{- 1} & 0 & 1 & 1 & 0 & 0 & 0 & 0 & 0 \\0 & {- 1} & {- 1} & 0 & 1 & 0 & 0 & 0 & 0 \\0 & 0 & 0 & {- 1} & {- 1} & 0 & 0 & 0 & 0 \\0 & 0 & 0 & 0 & 0 & 0 & {- 1} & 1 & 0 \\0 & 0 & 0 & 0 & 0 & {- 1} & 0 & 0 & 1 \\0 & 0 & 0 & 0 & 0 & 0 & 0 & {- 1} & 0 \\0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 1\end{bmatrix}_{9*9}$

In step S104, a sub-matrix that matches the first incidence matrix issearched from the third incidence matrix, that is, a sub-matrix thatmatches the incidence matrix of a topology substructure is searched fromthe incidence matrix of a primary topology structure of the power grid.FIG. 8 shows a flow chart of searching an incidence matrix of thetopology substructure from the incidence matrix of a primary topologystructure of the electric power grid. Specifically, the flow starts fromStep S801. In step S802, the row vectors and column vectors of theincidence matrix of a topology substructure is searched from the rowvectors and column vectors of the incidence matrix of a primary topologystructure, respectively. Specifically, FIG. 9 shows the specific processfor implementing Step S802 in FIG. 8. In step S8021, the row vectors ofthe incidence matrix of a topology substructure are searched one by onefrom the row vectors of the incidence matrix of a primary topologystructure. In step S8022, the row vectors searched out are added intothe row vectors' list, and the row numbers of the row vectors within theincidence matrix of the primary topology structure is recorded. In stepS8023, it is judged whether the row vectors' list includes all rowvectors in the incidence matrix of the primary topology structure. Ifresult of the judgment is no, the flow ends; if result of the judgmentis yes, the flow proceeds to Step S8024, in which the column vectors ofthe incidence matrix of a topology substructure is searched one by onefrom the column vectors of the incidence matrix of a primary topologystructure. In Step S8025, the column vectors searched out are added intothe column vectors' list, and the column numbers of the column vectorswithin the incidence matrix of the primary topology structure isrecorded. In step S8026, it is judged whether the column vectors' listincludes all column vectors in the incidence matrix of the primarytopology structure. If result of the judgment is no, the flow ends; ifresult of the judgment is yes, the flow proceeds to Step S803 in FIG. 8.Those skilled in the art will appreciate that when searching the rowvectors and column vectors of the incidence matrix of a topologysubstructure from the row vectors and column vectors of the incidencematrix of a primary topology structure, the sequence is not invariable.It is also possible to search column vectors of the incidence matrix ofa topology substructure first, and then search the row vectors of theincidence matrix of a topology substructure after making sure that allthe column vectors of the incidence matrix of a topology substructureare searched out. In this way, objective of this invention can also beachieved. In Step S803, if all the row vectors and column vectors of theincidence matrix of a topology substructure are searched out from therow vectors and column vectors of the incidence matrix of a primarytopology structure, it is determined that possibly matching sub-matrixcontained in the incidence matrix of primary topology structure. Byreferring to the row vectors and column vectors of the incidence matrixof a topology substructure, the row vectors and column vectors searchedout are permutated and combined, possibly matching sub-matrix containedin the incidence matrix of primary topology structure can be determined.In Step S804, a sub-matrix that matches the incidence matrix of thetopology substructure is filtered out from the possibly matchingsub-matrix contained in the incidence matrix of primary topologystructure. FIG. 10 shows a specific process for implementing Step S804in FIG. 8. In Step S8041, the vector coordinates of the possiblymatching sub-matrix are mapped with the vector coordinates of theincidence matrix of the topology substructure. In step S8042, the vectorvalues of the vector coordinates of the possibly matching sub-matrix iscompared with the vector values of the vector coordinates of theincidence matrix of the topology substructure having mapping relations.In step S8043, it is judged whether each vector value in the possiblymatching sub-matrix is same as the mapping vector values. If result ofthe judgment is no, then the flow proceeds to Step S805 and ends. Ifresult of the judgment is yes, then the flow proceeds to Step S8044, inwhich the sub-matrix that matches the incidence matrix of the topologysubstructure is filtered out, and the flow ends.

According to the above embodiment of this invention, the incidencematrix A of a topology substructure is searched from the incidencematrix A″ of a primary topology structure. With the searching method,the row vectors V of the incidence matrix A of the topology substructureis searched out one by one from the row vectors V′ of the incidencematrix A″. That is, the buses and dummy points are searched from rowvectors V′ of the incidence matrix A″. Because both auxiliary buses andmain buses are buses, the row vectors searched out include auxiliary bussegment 1, auxiliary bus segment 2, dummy point 3, dummy point 4, dummypoint 5, main bus segment 2, main bus segment 1. The row vectorssearched out are added into the row vectors' list as indicated in Table1, and the row number of the row vectors searched out within theincidence matrix of the primary topology structure is recorded. At thismoment, it is judged that the row vectors list includes all row vectorsof the incidence matrix of the topology substructure. Thus, the columnvectors M of the incidence matrix of the topology substructure issearched out one by one from the column vectors M″ of the incidencematrix of the primary topology structure. That is, the disconnector,breakers and ground disconnector are searched from the column vectors M″of the incidence matrix of the primary topology structure. The columnvectors searched out include disconnector 1, disconnector 2, breaker 1,ground disconnector 1, ground disconnector 2, breaker 2, breaker 3,disconnector 3, and disconnector 4. The column vectors searched out areadded into the column vectors' list as indicated in Table 2, and thecolumn number of the column vectors searched out within the incidencematrix of the primary topology structure is recorded. At this moment, itis judged that the column vectors list includes all column vectors ofthe incidence matrix of the topology substructure.

TABLE 1 Row vectors' list Row Vectors searched out Row No. auxiliary bussegment 1 1 auxiliary bus segment 2 2 main bus segment 1 7 main bussegment 2 6 dummy point 3 3 dummy point 4 4 dummy point 5 5

TABLE 2 Column vectors' list Column Vectors searched out Column No.disconnector 1 1 disconnector 2 2 disconnector 3 7 disconnector 4 6breaker 1 3 breaker 2 4 breaker 3 5 ground disconnector 1 4 grounddisconnector 2 5

Based on row vectors and column vectors searched out, the possiblymatching sub-matrix present in the incidence matrix of the topologysubstructure is determined. Based on principles of permutation andcombination, column vectors have the following 6 possibilities: U1, U2,. . . U6; row vectors have the following 18 possibilities: W1, W2, . . .W18. Thereby, ultimately there could be 108 possibly matchingsub-matrixes obtained: B1, B2, . . . B108. For the sake of briefdescription, only 3 of the possibilities, B1, B2 and B3, are given asexamples. Others are omitted herein.

$\mspace{20mu}{{{Row}\mspace{14mu}{vector}\mspace{14mu} U\; 1} = \left\lceil \begin{matrix}{{Auxiliary}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 1} \\{{Auxiliary}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 2} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 3} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 4} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 5}\end{matrix} \right\rbrack_{5*1}}$$\mspace{20mu}{{{Row}\mspace{14mu}{vector}\mspace{14mu} U\; 2} = \left\lceil \begin{matrix}{{Main}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 1} \\{{Main}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 2} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 3} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 4} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 5}\end{matrix} \right\rbrack_{5*1}}$$\mspace{20mu}{{{Row}\mspace{14mu}{vector}\mspace{14mu} U\; 3} = \begin{bmatrix}\begin{matrix}\begin{matrix}\begin{matrix}{{Main}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 1} \\{{Auxiliary}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 1}\end{matrix} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 3}\end{matrix} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 4}\end{matrix} \\{{{Dummy}\mspace{14mu}{point}\mspace{14mu} 5}\;}\end{bmatrix}_{5*1}}$$\mspace{20mu}{{{Row}\mspace{14mu}{vector}\mspace{14mu} U\; 4} = \begin{bmatrix}{{Main}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 1} \\{{Auxiliary}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 2} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 3} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 4} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 5}\end{bmatrix}_{5*1}}$$\mspace{20mu}{{{Row}\mspace{14mu}{vector}\mspace{14mu} U\; 5} = \begin{bmatrix}{{Main}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 2} \\{{Auxiliary}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 1} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 3} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 4} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 5}\end{bmatrix}_{5*1}}$$\mspace{20mu}{{{Row}\mspace{14mu}{vector}\mspace{14mu} U\; 6} = \begin{bmatrix}{{Main}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 2} \\{{Auxiliary}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 1} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 3} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 4} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 5}\end{bmatrix}_{5*1}}$$\mspace{20mu}{{{Column}\mspace{14mu}{vector}\mspace{14mu} W\; 1} = \begin{bmatrix}{{Disconnector}\mspace{14mu} 1} \\{{Disconnector}\mspace{14mu} 2} \\{{Breaker}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}_{5*1}^{T}}$$\mspace{20mu}{{{Column}\mspace{14mu}{vector}\mspace{14mu} W\; 2} = \begin{bmatrix}{{Disconnector}\mspace{14mu} 2} \\{{Disconnector}\mspace{14mu} 3} \\{{Breaker}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}_{5*1}^{T}}$$\mspace{20mu}{{{Column}\mspace{14mu}{vector}\mspace{14mu} W\; 3} = \begin{bmatrix}{{Disconnector}\mspace{14mu} 3} \\{{Disconnector}\mspace{14mu} 4} \\{{Breaker}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}_{5*1}^{T}}$$\mspace{20mu}{{{Column}\mspace{14mu}{vector}\mspace{14mu} W\; 4} = \begin{bmatrix}{{Disconnector}\mspace{14mu} 1} \\{{Disconnector}\mspace{14mu} 4} \\{{Breaker}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}_{5*1}^{T}}$$\mspace{20mu}{{{Column}\mspace{14mu}{vector}\mspace{14mu} W\; 5} = \begin{bmatrix}{{Disconnector}\mspace{14mu} 1} \\{{Disconnector}\mspace{14mu} 3} \\{{Breaker}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}_{5*1}^{T}}$$\mspace{20mu}{{{Column}\mspace{14mu}{vector}\mspace{14mu} W\; 6} = \begin{bmatrix}{{Disconnector}\mspace{14mu} 2} \\{{Disconnector}\mspace{14mu} 4} \\{{Breaker}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}_{5*1}^{T}}$$\mspace{20mu}{{{Column}\mspace{14mu}{vector}\mspace{14mu} W\; 7} = \begin{bmatrix}{{Disconnector}\mspace{14mu} 1} \\{{Disconnector}\mspace{14mu} 3} \\{{Breaker}\mspace{14mu} 2} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}_{5*1}^{T}}$$\mspace{20mu}{{{Column}\mspace{14mu}{vector}\mspace{14mu} W\; 8} = \begin{bmatrix}{{Disconnector}\mspace{14mu} 1} \\{{Disconnector}\mspace{14mu} 4} \\{{Breaker}\mspace{14mu} 2} \\{\;{{Ground}\mspace{11mu}{disconnector}\mspace{14mu} 1}} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}_{5*1}^{T}}$$\mspace{20mu}{{{Column}\mspace{14mu}{vector}\mspace{14mu} W\; 9} = \begin{bmatrix}{{Disconnector}\mspace{14mu} 2} \\{{Disconnector}\mspace{14mu} 3} \\{{Breaker}\mspace{14mu} 2} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}_{5*1}^{T}}$$\mspace{20mu}{{{Column}\mspace{14mu}{vector}\mspace{14mu} W\; 10} = \begin{bmatrix}{{Disconnector}\mspace{14mu} 2} \\{{Disconnector}\mspace{14mu} 4} \\{{Breaker}\mspace{14mu} 2} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}_{5*1}^{T}}$$\mspace{20mu}{{{Column}\mspace{14mu}{vector}\mspace{14mu} W\; 11} = \begin{bmatrix}{{Disconnector}\mspace{14mu} 3} \\{{Disconnector}\mspace{14mu} 4} \\{{Breaker}\mspace{14mu} 2} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}_{5*1}^{T}}$$\mspace{20mu}{{{Column}\mspace{14mu}{vector}\mspace{14mu} W\; 12} = \begin{bmatrix}{{Disconnector}\mspace{14mu} 1} \\{{Disconnector}\mspace{14mu} 2} \\{{Breaker}\mspace{14mu} 2} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}_{5*1}^{T}}$$\mspace{20mu}{{{Column}\mspace{14mu}{vector}\mspace{14mu} W\; 13} = \begin{bmatrix}{{Disconnector}\mspace{14mu} 1} \\{{Disconnector}\mspace{14mu} 3} \\{{Breaker}\mspace{14mu} 3} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}_{5*1}^{T}}$$\mspace{20mu}{{{Column}\mspace{14mu}{vector}\mspace{14mu} W\; 14} = \begin{bmatrix}{{Disconnector}\mspace{14mu} 1} \\{{Disconnector}\mspace{14mu} 2} \\{{Breaker}\mspace{14mu} 3} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}_{5*1}^{T}}$$\mspace{20mu}{{{Column}\mspace{14mu}{vector}\mspace{14mu} W\; 15} = \begin{bmatrix}{{Disconnector}\mspace{14mu} 1} \\{{Disconnector}\mspace{14mu} 4} \\{{Breaker}\mspace{14mu} 3} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}_{5*1}^{T}}$$\mspace{20mu}{{{Column}\mspace{14mu}{vector}\mspace{14mu} W\; 16} = \begin{bmatrix}{{Disconnector}\mspace{14mu} 2} \\{{Disconnector}\mspace{14mu} 3} \\{{Breaker}\mspace{14mu} 3} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}_{5*1}^{T}}$$\mspace{20mu}{{{Column}\mspace{14mu}{vector}\mspace{14mu} W\; 17} = \begin{bmatrix}{{Disconnector}\mspace{14mu} 3} \\{{Disconnector}\mspace{14mu} 4} \\{{Breaker}\mspace{14mu} 3} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}_{5*1}^{T}}$$\mspace{20mu}{{{Column}\mspace{14mu}{vector}\mspace{14mu} W\; 18} = \begin{bmatrix}{{Disconnector}\mspace{14mu} 2} \\{{Disconnector}\mspace{14mu} 4} \\{{Breaker}\mspace{14mu} 3} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}_{5*1}^{T}}$${B\; 1} = {{U\; 1 \times W\; 1} = {\begin{bmatrix}{{Auxiliary}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 1} \\{{Auxiliary}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 2} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 3} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 4} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 5}\end{bmatrix}_{5*1}\begin{bmatrix}{{Disconnector}\mspace{14mu} 1} \\{{Disconnector}\mspace{14mu} 2} \\{{Breaker}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}}_{5*1}^{T}}$${B\; 2} = {{U\; 1 \times W\; 2} = {\begin{bmatrix}{{Auxiliary}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 1} \\{{Auxiliary}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 2} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 3} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 4} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 5}\end{bmatrix}_{5*1}\begin{bmatrix}{{Disconnector}\mspace{14mu} 2} \\{{Disconnector}\mspace{14mu} 3} \\{{Breaker}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}}_{5*1}^{T}}$ ⋯${B\; 72} = {{U\; 6 \times W\; 18} = {\begin{bmatrix}{{Main}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 1} \\{{Auxiliary}\mspace{14mu}{bus}\mspace{14mu}{segment}\mspace{14mu} 2} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 3} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 4} \\{{Dummy}\mspace{14mu}{point}\mspace{14mu} 5}\end{bmatrix}_{5*1}\begin{bmatrix}{{Disconnector}\mspace{14mu} 2} \\{{Disconnector}\mspace{14mu} 4} \\{{Breaker}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 1} \\{{Ground}\mspace{14mu}{disconnector}\mspace{14mu} 2}\end{bmatrix}}_{5*1}^{T}}$

Tanking B1 as an example, according to method of generating incidencematrix, a sub-matrix B1 is obtained.

${B\; 1} = \begin{bmatrix}1 & 0 & 0 & 0 & 0 \\0 & 1 & 0 & 0 & 0 \\{- 1} & 0 & 1 & 1 & 0 \\0 & {- 1} & {- 1} & 0 & 1 \\0 & 0 & 0 & {- 1} & {- 1}\end{bmatrix}_{5*5}$

The vector coordinates of the possibly matching sub-matrix are mappedwith the vector coordinates of the incidence matrix of the topologysubstructure. Taking B1 as an example, the coordinates of vector b_(ij)of B1 are mapped with the coordinates of vector a_(ij) of the incidencematrix A of the topology substructure. By reference to the row vectors'table and column vectors' table, their mapping relations are recorded.

-   -   (b₁₁, a₁₁) (b₁₂, a₁₂) (b₁₃, a₁₃)(b₁₄, a₁₄)(b₁₅, a₁₅)    -   (b₂₁, a₂₁) (b₂₂, a₂₂) (b₂₃, a₂₃)(b₂₄, a₂₄)(b₂₅, a₂₅)    -   (b₃₁, a₃₁) (b₃₂, a₃₂) (b₃₃, a₃₃)(b₃₄, a₃₄)(b₃₅, a₃₅)    -   (b₄₁, a₄₁) (b₄₂, a₄₂) (b₄₃, a₄₃)(b₄₄, a₄₄)(b₄₅, a₄₅)    -   (b₅₁, a₅₁) (b₅₂, a₅₂) (b₅₃, a₅₃)(b₅₄, a₅₄)(b₅₅, a₅₅)

The vector values of vector coordinates having mapping relations, suchas the following values having mapping relations: (b₁₁, a₁₁), b₁₁=1,a₁₁=1, b₁₁=a₁₁; (b₁₂, a₁₂), b₁₂=0, a₁₂=0, b₁₂=a₁₂; . . . (b₅₅, a₅₅),b₅₅=−1, a₅₅=−1, b₅₅=a₅₅ are compared one by one. It is determinedtherefrom that each vector value of the sub-matrix B1 is same as thevector value of the mapped vector value in the incidence matrix A. Thenit is determined that the sub-structure B1 is a matching topologysubstructure.

Also taking B2 as an example, a substructure B2 is obtained by themethod of generating incidence matrix described above.

${B\; 2} = \begin{bmatrix}0 & 0 & 0 & 0 & 0 \\1 & 0 & 0 & 0 & 0 \\0 & 0 & 1 & 1 & 0 \\{- 1} & 0 & {- 1} & 0 & 1 \\0 & 0 & 0 & {- 1} & {- 1}\end{bmatrix}_{5*5}$

The coordinates of vector b_(ij) of B2 are mapped with the coordinatesof vector a_(ij) of the incidence matrix A of the topology substructure.By reference to the row vectors' table and column vectors' table, theirmapping relations are recorded.

-   -   (b₁₂, a₁₁) (b₁₈, a₁₂) (b₁₃, a₁₃)(b₁₄, a₁₄)(b₁₅, a₁₅)    -   (b₂₂, a₂₁) (b₂₈, a₂₂) (b₂₃, a₂₃)(b₂₄, a₂₄)(b₂₅, a₂₅)    -   (b₃₂, a₃₁) (b₃₈, a₃₂) (b₃₃, a₃₃)(b₃₄, a₃₄)(b₃₅, a₃₅)    -   (b₄₂, a₄₁) (b₄₈, a₄₂) (b₄₃, a₄₃)(b₄₄, a₄₄)(b₄₅, a₄₅)    -   (b₅₂, a₅₁) (b₅₈, a₅₂) (b₅₃, a₅₃)(b₅₄, a₅₄)(b₅₅, a₅₅)

The vector values of vector coordinates having mapping relations arecompared one by one. For example, if the following values having mappingrelations (b₁₂, a₁₁), b₁₂=0, a₁₁=1, b₁₂ and a₁₁ are not equal, it isdetermined therefrom that sub-matrix B2 is not a matching topologysubstructure.

Then, the remaining possibly matching sub-matrixes are traversed. It isfinally filtered out that only B1 is a matching topology substructure.The section surrounded by a black frame in the following incidencematrix A″ is the matching sub-matrix B1.

${{Incidence}\mspace{14mu}{Matrix}\mspace{14mu} A^{''}} = \left\lbrack {\begin{matrix}\begin{matrix}1 & 0 & 0 & 0 & 0 \\0 & 1 & 0 & 0 & 0 \\{- 1} & 0 & 1 & 1 & 0 \\0 & {- 1} & {- 1} & 0 & 1 \\0 & 0 & 0 & {- 1} & {- 1}\end{matrix} \\\begin{matrix}0 & \; & 0 & 0 & \; & 0 & 0 \\0 & \; & 0 & 0 & \; & 0 & 0 \\0 & \; & 0 & 0 & \; & 0 & 0 \\0 & \; & 0 & 0 & \; & 0 & 0\end{matrix}\end{matrix}\begin{matrix}\begin{matrix}1 & 0 & 0 & 0 \\0 & 1 & 0 & 0 \\0 & 0 & 0 & 0 \\0 & 0 & 0 & 0 \\0 & 0 & 0 & 0\end{matrix} \\\begin{matrix}0 & {- 1} & 1 & 0 \\{- 1} & 0 & 0 & 1 \\0 & 0 & {- 1} & 0 \\0 & 0 & 0 & {- 1}\end{matrix}\end{matrix}} \right\rbrack_{9*9}$

In step S105, a fourth incidence matrix is obtained by using the secondincidence matrix to transform the matching sub-matrix, in which thefourth incidence matrix corresponds to the incidence matrix of theskeleton topology structure of the primary topology structure.

In step S106, a skeleton topology structure corresponding to the primarytopology structure is generated according to the fourth incidencematrix.

According to an embodiment of this invention, obtaining the fourthincidence matrix by using the second incidence matrix to transform thematching sub-matrix comprises: replacing the matching sub-matrix withthe second incidence matrix to obtain the fourth incidence matrix; thengenerating the skeleton topology structure corresponding to the primarytopology structure according to the fourth incidence matrix.

According to the above embodiment of this invention, a incidence matrixB of skeleton topology structure corresponding to the primary topologystructure is generated by replacing the matching sub-matrix B1 with theincidence matrix A′ of the skeleton topology structure. The skeletontopology structure as shown in FIG. 11 corresponding to the primarytopology structure is generated according to the incidence matrix B ofthe skeleton topology structure corresponding to the primary topologystructure.

${{Incidence}\mspace{14mu}{Matrix}\mspace{14mu} A^{''}} = \left\lbrack {\begin{matrix}\begin{matrix}0 & 0 & {- 1} & 0 & 0 \\0 & 0 & 0 & 0 & 0 \\0 & 0 & 0 & 0 & 0 \\0 & 0 & 0 & 0 & 0 \\0 & 0 & 0 & {- 1} & {- 1}\end{matrix} \\\begin{matrix}0 & \; & 0 & 0 & \; & 0 & 0 \\0 & \; & 0 & 0 & \; & 0 & 0 \\0 & \; & 0 & 0 & \; & 0 & 0 \\0 & \; & 0 & 0 & \; & 0 & 0\end{matrix}\end{matrix}\begin{matrix}\begin{matrix}1 & 0 & 0 & 0 \\0 & 1 & 0 & 0 \\0 & 0 & 0 & 0 \\0 & 0 & 0 & 0 \\0 & 0 & 0 & 0\end{matrix} \\\begin{matrix}0 & {- 1} & 1 & 0 \\{- 1} & 0 & 0 & 1 \\0 & 0 & {- 1} & 0 \\0 & 0 & 0 & {- 1}\end{matrix}\end{matrix}} \right\rbrack_{9*9}$

According to another embodiment of this invention, generating a fourthincidence matrix by using the second incidence matrix to transform thematching sub-matrix comprises the steps of: determining a transformationmatrix according to the first incidence matrix and second incidencematrix; obtaining a fourth incidence matrix by using the transformationmatrix to transform the matching sub-matrix; generating a skeletontopology structure corresponding to the primary topology structureaccording to the fourth incidence matrix.

According to the above embodiment of this invention, taking incidencematrixes A and A′ computed in Step S102 as an example, thetransformation matrix is determined according to the incidence matrix ofthe topology substructure and the incidence matrix of the skeletontopology substructure. The relation among the incidence matrix of thetopology substructure, the incidence matrix of the skeleton topologysubstructure and the transformation matrix of the skeleton topologysubstructure is A×H=A′. Based on knowledge of matrix conversion, atransformation matrix is obtained with H=A⁻¹×A′. By using theformulation H=A⁻¹×A′, the transformation matrix H is obtained. Theincidence matrix corresponding to primary topology structure andskeleton topology structure is generated by using the transformationmatrix H and matching sub-matrix B1. By using formulation B1×H=A′, B1 istransformed into A′, and incidence matrix B of skeleton topologystructure corresponding to the primary topology structure is generated.

$\mspace{20mu}{{{Transformation}\mspace{14mu}{matrix}\mspace{14mu} H} = \begin{bmatrix}0 & 0 & 1 & 0 & 0 \\0 & 0 & {- 1} & 0 & 0 \\0 & 0 & 1 & 0 & 0 \\0 & 0 & 0 & 0 & 0 \\0 & 0 & 0 & 0 & 0\end{bmatrix}_{5*5}}$$\mspace{20mu}{A^{\prime} = {{B\; 1 \times H} = {\begin{bmatrix}1 & 0 & 0 & 0 & 0 \\0 & 1 & 0 & 0 & 0 \\{- 1} & 0 & 1 & 1 & 0 \\0 & {- 1} & {- 1} & 0 & 1 \\0 & 0 & 0 & {- 1} & {- 1}\end{bmatrix}_{5*5}\begin{bmatrix}0 & 0 & 1 & 0 & 0 \\0 & 0 & {- 1} & 0 & 0 \\0 & 0 & 1 & 0 & 0 \\0 & 0 & 0 & 0 & 0 \\0 & 0 & 0 & 0 & 0\end{bmatrix}}_{5*5}}}$${{Transformation}\mspace{14mu}{m{atrix}}\mspace{14mu} A^{''}} = \left\lbrack {\begin{matrix}\begin{matrix}1 & 0 & 0 & 0 & 0 \\0 & 1 & 0 & 0 & 0 \\{- 1} & 0 & 1 & 1 & 0 \\0 & {- 1} & {- 1} & 0 & 1 \\0 & 0 & 0 & {- 1} & {- 1}\end{matrix} \\\begin{matrix}0 & \; & 0 & 0 & \; & 0 & 0 \\0 & \; & 0 & 0 & \; & 0 & 0 \\0 & \; & 0 & 0 & \; & 0 & 0 \\0 & \; & 0 & 0 & \; & 0 & 0\end{matrix}\end{matrix}\begin{matrix}\begin{matrix}1 & 0 & 0 & 0 \\0 & 1 & 0 & 0 \\0 & 0 & 0 & 0 \\0 & 0 & 0 & 0 \\0 & 0 & 0 & 0\end{matrix} \\\begin{matrix}0 & {- 1} & 1 & 0 \\{- 1} & 0 & 0 & 1 \\0 & 0 & {- 1} & 0 \\0 & 0 & 0 & {- 1}\end{matrix}\end{matrix}} \right\rbrack_{9*9}$

Based on the same inventive concept, this invention provides a devicefor extracting skeleton topology of electric power grid. FIG. 12 shows ablock diagram of the device for extracting skeleton topology in anelectric power grid, according to an embodiment of this invention. Thedevice comprises: a receiving module 1201 configured to receive adescription of a topology sub-structure corresponding with user's needand a description of skeleton topology sub-structure extracted from thetopology sub-structure; an incidence matrix generating module 1202configured to generate a first incidence matrix based on the descriptionof the topology sub-structure, to generate a second incidence matrixbased on the description of the skeleton topology sub-structure, and togenerate a third incidence matrix based on a primary topology structureof the electric power grid; a searching module 1203 configured to searchfrom the third incidence matrix a sub-matrix that matches the firstincidence matrix; a matrix transforming module 1204 configured to obtaina fourth incidence matrix by using the second incidence matrix totransform the matching sub-matrix; and a skeleton topology structuregenerating module 1205 configured to generate the skeleton topologystructure corresponding to the primary topology structure according tothe fourth incidence matrix.

According to an embodiment of this invention, the matrix transformingmodule 1204 is further configured to: generate a fourth incidence matrixby using the second incidence matrix to replace the matching sub-matrix.

According to an embodiment of this invention, the matrix transformingmodule 1204 is further configured to: determine the transformationmatrix according to the first incidence matrix and second incidencematrix; and obtain a fourth incidence matrix by using the transformationmatrix to transform the matching sub-matrix.

According to an embodiment of this invention, the searching module 1203further comprises: a vector searching module, being configured to searchthe row vectors and column vectors of the first incidence matrix fromthe row vectors and column vectors of the third incidence matrixrespectively; a sub-matrix determining module, being configured todetermine that possibly matching sub-matrix is contained in the thirdincidence matrix if all the row vectors and column vectors of the firstincidence matrix are searched out from the row vectors and columnvectors of the third incidence matrix; a sub-matrix filtering module,being configured to filter out the sub-matrix matching the firstincidence matrix from the possibly matching sub-matrixes contained inthe third incidence matrix.

According to an embodiment of this invention, the vector searchingmodule is further configured to: search the row vectors of the firstincidence matrix one by one from the row vectors of the third incidencematrix; add the row vectors searched out into the row vectors' list, andrecord the row number of the row vectors in the third incidence matrix;search the column vectors of the first incidence matrix from the columnvectors of the third incidence matrix if it is determined that the rowvectors' list includes all row vectors in the first incidence matrix;add the column vectors searched out into the column vectors' list, andrecord the column number of the column vectors in the third incidencematrix; judge whether the column vectors' list includes all columnvectors in the first incidence matrix;

According to an embodiment of this invention, the sub-matrix determiningmodule is further configured to: by reference to the row vectors andcolumn vectors of the incidence matrix of a topology substructure,permutate and combine the row vectors and column vectors searched out todetermine possibly matching sub-matrix contained in the third incidencematrix.

According to an embodiment of this invention, the sub-matrix filteringmodule is further configured to: map the vector coordinates of thepossibly matching sub-matrix with the vector coordinates of the firstincidence matrix; to compare the vector values of the vector coordinatesof the possibly matching sub-matrix with the vector values of the vectorcoordinates of the first incidence matrix having mapping relations; ifthe vector value of each vector coordinate of the possibly matchingsub-matrix is same as the vector value of vector coordinates of thefirst incidence matrix having mapping relations, to filter out thesub-matrix matching the first incidence matrix from the possiblymatching sub-matrixes.

The method and device for extracting skeleton topology structureaccording to this invention can automatically and efficiently extractskeleton topology structures according to business needs. Thus, thefrequent occurrence of collisions caused by different frequencies inupdating skeleton topology structures by various business departmentscan be avoided.

FIG. 13 is a schematic block diagram of the structure of the computingdevices for implementing the embodiments of this invention. The computersystem in FIG. 13 comprises CPU (central processing unit) 1301, RAM(random access memory) 1302, ROM (read-only memory) 1303, system bus1304, hard disc controller 1305, keyboard controller 1306, serialinterface controller 1307, parallel interface controller 1308, displaycontroller 1309, hard disc 1310, keyboard 1311, serial peripheralequipment 1312, parallel peripheral equipment 1313 and display 1314.Among these components, the CPU 1301, RAM 1302, ROM 1303, hard disccontroller 1305, keyboard controller 1306, serial interface controller1307, parallel interface controller 1308, display controller 1309 areconnected to system bus 1304; hard disc 1310 is connected to hard disccontroller 1305; keyboard 1311 is connected to keyboard controller 1306;serial peripheral equipment 1312 is connected to serial interfacecontroller 1307; parallel peripheral equipment 1313 is connected toparallel interface controller 1308; and display 1314 is connected todisplay controller 1309.

The function of each component in FIG. 13 is well known in this field,and the structure in FIG. 13 is conventional. Such a structure is notonly used in PC, but also in handhold devices such as Palm PC, PDA(personal digital assistant), mobile phones. In different appliances,such as in realizing the user terminal containing the user moduleaccording to this invention or the host server containing the networkapplication servers according to this invention, some components can beadded to the structure shown in FIG. 13, or some components in FIG. 13can be omitted. The system shown in FIG. 13 is controlled by computerreadable instructions, generally stored in hard disc 1310, EPROM orother nonvolatile memories as software. The software can also bedownloaded from internet (not indicated in the drawing). The softwarestored in hard disc 1310 or downloaded from internet can be loaded inRAM 1302 and executed by CPU 1301, so as to carry out functions definedby the software.

Although the computer system described in FIG. 13 can support thetechnical solution provided according to this invention, it is just anexample of a computer system. Those skilled in the art will appreciatethat embodiments of this invention can also be realized with many otherdesigns of computer systems.

The exemplarily embodiments of the present invention have been describedby reference to the attached drawings, but it should be understood thatthis invention is not limited to these exact embodiments, and those ofordinary skill in the art can make modifications and variations withoutdeparting from the scope and spirit of the invention. All thesemodifications and variations are intended to be included within thescope of this invention defined by the attached claims.

It should be understood that at least some aspects of this invention canbe alternatively implemented as program products. Programs definingfunctions relating to this invention can be transmitted to data storagedevice or computer devices through various signal bearing media. Thesignal bearing media include but not limit to read only medium (such asCD-ROM), writable storage medium (such as floppy disc and hard discdriver, read/write CD ROM, optical medium) as well as communicationmedia like computers having ethernet and telephone networks. Thereby, itshould be understood that among these signal bearing media, when theycarry or coded with the computer readable instructions managing themethodology functions of this invention, they represent alternativeembodiments of this invention. This invention can be realized in formsof hardware, software, firmware or their combinations. This inventioncan be realized on one computer device in a centralized manner, or berealized in a dispersed manner. With the dispersed manner, differentcomponents are dispersed on a number of inter-connected computerdevices. Any computer devices or other devices that are suitable forexecuting the methods described herein are appropriate. Preferably, thisinvention is carried out in a form of combination of computer softwareand general purpose computer hardware. In this way of realizing theinvention, when the computer program is loaded and executed, thecomputer device is controlled to execute the method of this invention orconstitute the system of this invention.

The above description is intended to describe with examples thepreferred embodiment of this invention. The above descriptions ofpreferred embodiments are not exhaustive, neither it is intended tolimited to the invention in the definite form disclosed. In view of thisteaching, many modifications and variations are possible. Thesemodifications and variations are apparent for those skilled in the art,and are within the scope of this invention defined by the attachedclaims.

What is claimed is:
 1. A computer-implemented method for extracting askeleton topology for an electric power grid, comprising: providing, byone or more processors of a computer, a primary topology structure ofthe electric power grid to a user via a graphic interface of thecomputer; receiving, by the one or more processors via an editing toolof the graphic interface that is operable by the user, a selection of aselected topology substructure from the primary topology structure and askeleton topology substructure extracted from the selected topologysubstructure; generating, by the one or more processors, a firstincidence matrix based on the selected topology substructure, the firstincidence matrix indicating connections between nodes representingelectrical components within the selected topology substructure, thefirst incidence matrix indicating inflow and outflow of each connectionof each node; generating, by the one or more processors, a secondincidence matrix based on the skeleton topology substructure, the secondincidence matrix indicating connections between nodes representingelectrical components within the skeleton topology substructure, thesecond incidence matrix indicating inflow and outflow of each connectionof each node; generating, by the one or more processors, a thirdincidence matrix based on the primary topology structure, the thirdincidence matrix indicating connections between nodes representingelectrical components within the primary topology structure, the thirdincidence matrix indicating inflow and outflow of each connection ofeach node; searching, by the one or more processors, from the thirdincidence matrix a sub-matrix that matches the first incidence matrix,the sub-matrix indicating connections between nodes representingelectrical components within the selected topology substructure, thesub-matrix indicating inflow and outflow of each connection of eachnode; obtaining, by the one or more processors, a fourth incidencematrix by using the second incidence matrix to transform the sub-matrix,the fourth incidence matrix indicating connections between nodesrepresenting electrical components within a skeleton topology structure,the fourth incidence matrix indicating inflow and outflow of eachconnection of each node; and generating, by the one or more processors,the skeleton topology structure extracted from the selected topologysubstructure based on the fourth incidence matrix, the skeleton topologystructure enabling the user to control effectiveness of the electricalpower grid.
 2. The computer-implemented method according to claim 1,wherein obtaining the fourth incidence matrix by using the secondincidence matrix to transform the sub-matrix comprises: replacing, bythe one or more processors, the sub-matrix with the second incidencematrix to generate the fourth incidence matrix.
 3. Thecomputer-implemented method according to claim 1, wherein obtaining afourth incidence matrix by using the second incidence matrix totransform the sub-matrix comprises: determining, by the one or moreprocessors, a transformation matrix according to the first incidencematrix and the second incidence matrix; obtaining, by the one or moreprocessors, the fourth incidence matrix by using the transformationmatrix to transform the sub-matrix.
 4. The computer-implemented methodaccording to claim 1, wherein searching from the third incidence matrixa submatrix that matches the first incidence matrix comprises:searching, by the one or more processors, row vectors and column vectorsof the first incidence matrix from row vectors and column vectors of thethird incidence matrix respectively; when all row vectors and columnvectors of the first incidence matrix are searched out from the rowvectors and column vectors of the third incidence matrix, determining,by the one or more processors, a possibly matching submatrix containedin the third incidence matrix based on the row vectors and columnvectors searched out; and filtering out, by the one or more processors,the sub-matrix that matches the first incidence matrix from the possiblymatching sub-matrix contained in the third incidence matrix.
 5. Thecomputer-implemented method according to claim 4, wherein searching rowvectors and column vectors of the first incidence matrix from rowvectors and column vectors of the third incidence matrix respectivelycomprises: searching, by the one or more processors, the row vectors ofthe first incidence matrix one by one from the row vectors of the thirdincidence matrix; adding, by the one or more processors, the row vectorssearched out into a row vectors' list, and recording a row number foreach of the row vectors searched out in the third incidence matrix; whendetermining that the row vectors' list includes all row vectors in thefirst incidence matrix, searching, by the one or more processors, thecolumn vectors of the first incidence matrix from the column vectors ofthe third incidence matrix one by one; adding, by the one or moreprocessors, the column vectors searched out into a column vectors' list,and recording a column number for each of the column vectors in thethird incidence matrix; and judging, by the one or more processors,whether the column vectors' list includes all column vectors in thefirst incidence matrix.
 6. The computer-implemented method according toclaim 4, wherein determining the possibly matching sub-matrix containedin the third incidence matrix based on the row vectors and columnvectors searched out comprises: determining, by the one or moreprocessors, the possibly matching submatrix contained in the thirdincidence matrix by permutating and combining the row vectors and columnvectors searched out by reference to the row vectors and column vectorsof the first incidence matrix.
 7. The computer-implemented methodaccording to claim 6, wherein filtering out the sub-matrix matching thefirst incidence matrix from the possibly matching sub-matrix containedin the third incidence matrix comprises: mapping, by the one or moreprocessors, vector coordinates of the possibly matching sub-matrix withvector coordinates of the first incidence matrix; comparing, by the oneor more processors, vector values of the vector coordinates of thepossibly matching sub-matrix with vector values of the vectorcoordinates of the first incidence matrix having mapping relations; andwhen the vector value of each vector coordinate of the possibly matchingsub-matrix is same as the vector value of vector coordinates of thefirst incidence matrix having mapping relations, filtering out, by theone or more processors, the sub-matrix matching the first incidencematrix from the possibly matching sub-matrix.
 8. A device for extractinga skeleton topology for an electric power grid, comprising: one or moreprocessors; and a memory, coupled to the one or more processors, thememory storing processor-executable program instructions, wherein theone or more processors, when executing the program instructions, areconfigured to: provide a primary topology structure of the electricpower grid to a user via a graphic interface of the device; receive viaan editing tool of the graphic interface that is operable by the user aselection of a selected topology substructure from the primary topologystructure and a skeleton topology substructure extracted from theselected topology substructure; generate a first incidence matrix basedon the selected topology substructure, the first incidence matrixindicating connections between nodes representing electrical componentswithin the selected topology substructure, the first incidence matrixindicating inflow and outflow of each connection of each node; generatea second incidence matrix based on the skeleton topology substructure,the second incidence matrix indicating connections between nodesrepresenting electrical components within the skeleton topologysubstructure, the second incidence matrix indicating inflow and outflowof each connection of each node; generate a third incidence matrix basedon the primary topology structure, the third incidence matrix indicatingconnections between nodes representing electrical components within theprimary topology structure, the third incidence matrix indicating inflowand outflow of each connection of each node; search from the thirdincidence matrix a sub-matrix that matches the first incidence matrix,the sub-matrix indicating connections between nodes representingelectrical components within the selected topology substructure, thesub-matrix indicating inflow and output of each connection of each node;obtain a fourth incidence matrix by using the second incidence matrix totransform the sub-matrix, the fourth incidence matrix indicatingconnections between nodes representing electrical components within askeleton topology structure, the fourth incidence matrix indicatinginflow and outflow of each connection of each node; and generate theskeleton topology extracted from the selected topology substructurebased on the fourth incidence matrix, the skeleton topology enabling theuser to control effectiveness of the electrical power grid.
 9. Thedevice according to claim 8, wherein to obtain the fourth incidencematrix, the one or more processors are further configured to: generatethe fourth incidence matrix by using the second incidence matrix toreplace the sub-matrix.
 10. The device according to claim 8, wherein toobtain a fourth incidence matrix, the one or more processors are furtherconfigured to: determine a transformation matrix according to the firstincidence matrix and the second incidence matrix; and obtain the fourthincidence matrix by using the transformation matrix to transform thesub-matrix.
 11. The device according to claim 8, wherein to search fromthe third incidence matrix a sub-matrix, the one or more processors arefurther configured to: search row vectors and column vectors of thefirst incidence matrix from row vectors and column vectors of the thirdincidence matrix respectively; determine a possibly matching sub-matrixcontained in the third incidence matrix when all the row vectors andcolumn vectors of the first incidence matrix are searched out from therow vectors and column vectors of the third incidence matrix; and filterout the sub-matrix matching that matches the first incidence matrix fromthe possibly matching sub-matrix contained in the third incidencematrix.
 12. The device according to claim 11, wherein to search the rowvectors and column vectors of the first incidence matrix, the one ormore processors are further configured to: search the row vectors of thefirst incidence matrix one by one from the row vectors of the thirdincidence matrix; add the row vectors searched out into a row vectors'list, and record a row number for each of the row vectors searched outin the third incidence matrix; when determining that the row vectors'list includes all row vectors in the first incidence matrix, searchcolumn vectors of the first incidence matrix one by one from the columnvectors of the third incidence matrix; add the column vectors searchedout into a column vectors' list, and record a column number for each ofthe column vectors in the third incidence matrix; and judge whether thecolumn vectors' list includes all column vectors of the first incidencematrix.
 13. The device according to claim 11, wherein to determine apossibly matching sub-matrix contained in the third incidence matrix,the one or more processors are further configured to: determine thepossibly matching sub-matrix contained in the third incidence matrix bypermutating and combining the row vectors and column vectors searchedout by reference to the row vectors and column vectors of the firstincidence matrix.
 14. The device according to claim 13, wherein tofilter out the sub-matrix matching the first incidence matrix from thepossibly matching sub-matrix, the one or more processors are furtherconfigured to: map vector coordinates of the possibly matchingsub-matrix with vector coordinates of the first incidence matrix;compare vector values of the vector coordinates of the possibly matchingsub-matrix with vector values of the vector coordinates of the firstincidence matrix having mapping relations; and when the vector value oreach vector coordinate of the possibly matching sub-matrix is same asthe vector value of vector coordinates of the first incidence matrixhaving mapping relations, filter out the sub-matrix matching the firstincidence matrix from the possibly matching submatrix.
 15. A computerprogram product comprising a non-transitory computer-readable storagemedium having computer readable program instructions embodied therewith,the computer readable program instructions executed by one or moreprocessors to cause a computer to perform a method for extracting askeleton topology for an electric power grid, said method comprising:providing, by the one or more processors of the computer, a primarytopology structure of the electric power grid to a user via a graphicinterface of the computer; receiving, by the one or more processors viaan editing tool of the graphic interface that is operable by the user, aselection of a selected topology substructure from the primary topologystructure and a skeleton topology substructure extracted from theselected topology substructure; generating, by the one or moreprocessors, a first incidence matrix based on the topology substructure,the first incidence matrix indicating connections between nodesrepresenting electrical components within the topology substructure, thefirst incidence matrix indicating inflow and outflow of each connectionof each node; generating, by the one or more processors, a secondincidence matrix based on the skeleton topology substructure, the secondincidence matrix indicating connections between nodes representingelectrical components within the skeleton topology substructure, thesecond incidence matrix indicting inflow and outflow of each connectionof each node; generating, by the one or more processors, a thirdincidence matrix based on the primary topology structure, the thirdincidence matrix indicating connections between nodes representingelectrical components within the primary topology structure, the thirdincidence matrix indicating inflow and outflow of each connection ofeach node; searching, by the one or more processors, from the thirdincidence matrix a sub-matrix that matches the first incidence matrix,the sub-matrix indicating connections between nodes representingelectrical components within the selected topology substructure, thesub-matrix indicating inflow and outflow of each connection of eachnode; obtaining, by the one or more processors, a fourth incidencematrix by using the second incidence matrix to transform the sub-matrix,the fourth incidence matrix indicating connections between nodesrepresenting electrical components within a skeleton topology structure,the fourth incidence matrix indicating inflow and outflow of eachconnection of each node; and generating, by the one or more processors,the skeleton topology structure extracted from the selected topologysubstructure based on the fourth incidence matrix, the skeleton topologystructure enabling the user to control effectiveness of the electricalpower grid.